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Abstract 

We give a new proof for a Ray-Knight representation of Feller's 
branching diffusion with logistic growth in terms of the local times 
of a reflected Brownian motion H with a drift that is afflne linear 
in the local time accumulated by H at its current level. In [5], such 
a representation was obtained by an approximation through Harris 
paths that code the genealogies of particle systems. The present proof 
is purely in terms of stochastic analysis, and is inspired by previous 
work of Norris, Rogers and Williams [7]. 

1 Introduction 

The second one of the two classical Ray-Knight theorems (see e.g. [TU] or 
establishes a representation of Feller's branching diffusion in terms of 
reflected Brownian motion. In words, it may be stated as follows: Take 
a standard Brownian motion on M + reflected at 0, and stopped when the 
local time accumulated at reaches a value x. Then the (total) local time 
accumulated by the resulting path at "height" t, viewed as a process indexed 
by t, is a Feller branching diffusion obeying the SDE dZf = 2y^ZfdWf 
with Zq = x. One way to interpret this is to view the reflected Brownian 
path as an exploration path which codes the genealogy of a continuous state 
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branching process (see e.g. [6]): the local time of the exploration path at 
height t measures the "width of the genealogical forest" at this level, or 
equivalently, the mass (or size) of the population that is alive at the time 
corresponding to this height. This mass is Zt, the state of the branching 
process at time t. The exploration path is a concatenation of Brownian 
excursions, with each excursion corresponding to a continuum random tree 
in the sense of Aldous [I]. The independence of the excursions leads to 
independent increments of (Z x ) when viewed as a process indexed by x. 
This mutual independence of offspring coming from different ancestors is 
also referred to as the branching property. 

Consider now, instead of a Feller branching diffusion, the same SDE with 
a logistic drift, namely 



where #,7 > are given parameters which will be fixed throughout the 
paper, and x > is the initial population. Just like Feller's branching 
diffusion, also the solution of (pQ), called Feller's branching diffusion with 
logistic growth, arises as the diffusion limit of discrete population models, 
but now with an interaction between individuals. This interaction can be 
thought of as a competition among individuals for resources, resulting in 
"pairwise lethal fights" (with intensity 7) that counteract the supercritical 
growth of the population. This population model and its diffusion limit ([T|) 
have been studied in j3j. 

In this note we will specify an SDE for a process (H s ), from which Feller's 
branching diffusion with logistic growth can be read off in the same way as 
Feller's branching diffusion is read off from reflected Brownian motion. Our 
proof will rely purely on stochastic analysis. Still, we give here some brief 
explanations on the underlying population model (for more illustrations and 
background we refer to our survey paper [8j). 

The process (H s ) will be reflected Brownian motion with a drift that 
depends on the local time t accumulated at the current level H s up to time 
s. More specifically, the drift coefficient will be of the form 9/2— jl. One way 
to understand the form of the drift is to see (H s ) again as the exploration 
process of a forest of random real trees, and to think of an approximation in 
terms of piecewise linear, continuous processes with constant absolute slope 
(so-called Harris paths): the rate of minima (giving rise to new branches) is 
increased proportional to 9/2, and the rate of maxima (describing deaths of 
branches) is increased proportional to the number of individuals visited by 
the exploration process so far on the current level. In our recent work |5j 
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we obtained the process H as the limit of exploration processes of discrete 
population models, and in this way provided a Ray-Knight representation 
of dU). The discrete approach of [5] gives a worthy insight into the result, 
since the way how the genealogy is built and how the exploration process 
codes the genealogical tree of the population is readily understandable at 
the discrete level. 

The derivation presented in this note does not rely on a discrete ap- 
proximation, but directly exploits methods from stochastic analysis. We 
use ideas from previous work of Norris, Rogers and Williams [7J to which 
our attention was drawn after the completion of [5] thanks to a hint of J-F 
Le Gall. In [7] a generalization of the first Ray-Knight theorem for "Brown- 
ian motions with a local time drift" was provided for cases that include the 
drift appearing in the SDE ([5]). 

We also extend the Ray-Knight representation of ([1]) by establishing an 
equality between laws of random fields (random functions of time t and 
ancestral mass x). In Section [2] we introduce Feller's branching diffusion 
with logistic growth as a random field {Zf, t,x > 0}. This is a natural set- 
up for the formulation of our main result, which is given in Section [3] and 
whose proof is contained in Section HJ The last section gives two remarks 
concerning a possible shortcut in the proof of the Theorem, and a general 
version of the second Ray-Knight theorem in the framework of [7] . 

2 A coupling over the ancestral masses 

In this section we define a random field {Zf, t,x > 0} such that for any 
x > 0, Z x := {Zf, t > 0} is a Feller branching diffusion with logistic 
growth and ancestral mass x, for any t > 0, x i— ¥ Zf is non-decreasing and 
x i — y {Zf, t > 0} is a (path- valued) Markov process which will be specified 
below.. 

To this purpose we define a family of transition probabilities Px,x > 0, 
on E, where E = C C (M+,M+) is the set of continuous mappings from M + 
to M + with compact support. Here and below, R + = [0, +oo). For x > 
and z e C C (IR + ,1R + ), let P x (z, •) be the distribution of z + Z z > x , where Z z > x 
solves 

Z x t ' z = x + f ZY{0 - 7 [Z^ Z + 2z(u)])du + 2 f ^/z^dWu, (2) 
Jo Jo 

with W being a standard Brownian motion. The equality P x (z,E) = 1 
is valid because Z x, ° (and a fortiori Z x,z ) a.s. hits zero in finite time (for 
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a proof of this fact see e.g. [3]). Before we show in Lemma Q] that the 
transition probabilities P x ,x > 0, indeed specify a random field {Zf}, let 
us briefly motivate the form of (|2|) in the light of the interpretation given in 
the introduction. 

Consider the evolution of the progeny of a sum of ancestral masses z(0) + 
x. The offspring of z(0) is assumed to follow (in an autonomous way) the 
dynamics of a Feller branching with logistic drift, the path z(t) stands for a 
realization of this offspring. The progeny (Z^' z ) of the additional ancestral 
mass x does not evolve independently of the offspring of x, but experiences 
an additional pressure coming from the given z(t), resulting in the negative 
drift — 2^z{t)Z^' z . In a finite population approximation with k+i ancestors, 
this means that the descendants of the £ "additional" ancestors suffer from 
the competition with those of the "first" k, while the descendants of the 
"first" k do not feel the presence of the descendants of the I "additional" 
ancestors. Recall from the introduction that in an individual-based model 
the competition leading to the negative quadratic drift in the populaiton 
size is modelled by pairwise fights between the individuals. If we think of 
the individuals being arranged in a linear order "from left to right" , where 
this order is passed on to the individual's offspring, then in our convention 
the pairwise fights are always be won by the individual to the left, resulting 
in the death of the individual to the right. 

Lemma 1 The family P x ,% > 0, satisfies the Chapman- Kolmogorov rela- 
tions. 

Proof: Observe that conditioned on Z x ' z , the random path V := Z V)Z+ZX ' Z 



with W being a standard Brownian motion (independent of W). Note that 
in the case 7 = 0, the two processes Z x,z and V are independent, as they 
should. Now Z x ' z + V satisfies 



This shows that z + Z x,z + V has distribution P x+y (z, •), as required. Indeed, 
since the two Brownian motions W and W' are independent, 



solves 




(3) 
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and consequently, from a well-known martingale representation theorem, 
there exists a third standard Brownain motion W" such that 




Definition 1 Let {Z x } x >o be the C C (M + ,1R + ) -valued Markov process with 
transition semigroup (P x ), starting from Z° = 0, the null trajectory. 

Remark 1 For each x > 0, Z x solves the SDE 

dZ' x = [9Zf - j(Z? ) 2 ] dt + 2yfZ*dW?, Z* = x, (4) 

where {W®, t > 0} is a standard Brownian motion. Since for x,y > the 
increment Z x+y — Z x is driven by a Brownian motion independent of that 
driving Z x , we have d(Z x , Z x+y ) t = d(Z x ,Z x ) t = Zf dt and conseqently 

d{W x , W x+y ) t = ^Z x /Zt +V dt, with the convention § = 0. 

3 A Ray-Knight representation 

Consider the following SDE driven by standard Brownian motion B 

H S = B S + \l s {Q) + e -s-~i f L r {H r )dr, s > 0, (5) 
1 1 Jo 

Here and everywhere below, {L s (t), s > 0, t > 0} denotes the local time of 
the process {H s , s > 0} accumulated up to time s at level t. Proposition [21 
stated and proved in the next section, will ensure (by specializing it to the 
case z = 0) that equation ([5]) has a unique weak solution, which we assume 
to be defined on some probability space (Q, J 7 , P). 
Define for any x > the stopping time 

S x = inf{s > 0, L s (0) > x}, 

and let {Zf, x,t > 0} denote the random field constructed in Section [2j 

Our main result is the 

Theorem The two random fields {Ls x (t), t, x > 0} and {Zf, t,x > 0} 
have the same law. 
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4 Proof of the Theorem 

To prepare for the proof of the Theorem, we first fix a z S C C (M+,IR+) and 
consider the SDE 

R* a =B S + l -L z M + °-s -7 J\< H r) + WH?)}dr, s > 0, (6) 

where L z stands for the local time of H z . We will prove in Subsection 14.11 

Proposition 2 The SDE (J5J) has a unique weak solution. 

Suppressing the superscript z, define for any x > the stopping time 

S x = inf{s > 0, L*(0) > x}. (7) 

The main step in the proof of the Theorem will be to show 

Proposition 3 For x > and z 6 C C (M+,IR + ) let {Z^ z , t > 0} be the 
solution of ([2|). Then the two processes {L z s (t), t > 0} and {Z^' z , t > 0} 
have the same law. 

4.1 Proof of Proposition [2] 

Let H denote Brownian motion reflected above 0, i. e. 

H S = B S + ~L S (0), 

where B is a J-^-standard Brownian motion defined on a probability space 
(O, J 7 , P), with J 7 = J 7 ^, and L is the semimartingale local time of H. Let 

G s = exp(M s -^(M) s Y s>0, 

with M s := J S {f - 7 [z(IT r ) + L r (i? r )]} dB r . (Recall that z G C C (M + ,M + ) 
is fixed.) The condition 

E(G S ) = 1, Vs > (8) 

is sufficient for the local martingale {G s , s > 0} to be a martingale. In that 
case, there exists a new probability measure P on such that for all 

s > 0, 
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and it follows from Girsanov's theorem (see e. g. Theorem VIII 1.4 in |10j ) 
that 

B a :=B a - J !^-j[z(H r ) + L r (H r )]\dr, s > 0, (9) 

is a standard Brownian motion. (Note that this does not require that P 
be absolute continuity with respect to P on J 7 .) Hence existence of a weak 
solution to © follows from 0, which in turn (see Theorem 1.1, chapter 7, 
page 152, in [3j) follows if we can ensure that for each s > there exists 
constants a > such that 

sup Eexp(ai? r ) < oo, (10) 

0<r<s 

where R r = || — j[z(H r ) + L r (H r )]\ 2 . Since z is bounded, the inequality 



pop is immediate from the following 

Lemma 2 Let H be a Brownian motion on M + reflected at the origin. Then 
for all s > there exists a = a(s) > such that 

sup E (exp(aL r (ff r ) 2 )) < oo. 

0<r<s 

Proof: Together with a simple scaling argument and a desintegration with 
respect to H r , this is immediate from the following 

Lemma 3 Let /3 be a standard Brownian motion starting at 0, and denote 
by Li(y) the local time accumulated by \f3\ at position y up to time 1. There 
exist constants a > and c > (not depending on y) such that for almost 
ally>0 

E[ e aL i(f) 2 | | ft | = y ] < c. (11) 

Proof: By symmetry the l.h.s. of ([TT]) a.s. equals E[e aLl ^ 2 | ft = y]. 
Writing ¥ y for the probability measure of a Brownian bridge from the origin 
at time to position y at time 1, and K\{a) for the local time accumulated 
up to time 1 at position a, we thus have to show the inequality 

for suitable constants a and c. By the Cauchy-Schwarz inequality, the l.h.s. 
of (fT2j) is bounded by 



1/2 



E- y 



1/2 



(13) 
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To estimate the first factor, we desintegrate with respect to the time U at 
which the path of the Brownian bridge first hits the level y. Conditioned 
under {U = u}, the part before time u does not contribute to the local time 
at y, and the second part is (by the stong Markov property) a Brownian 
bridge from y to y over a time interval of length 1 — u, hence (by scaling) 
the distribution of its local time at y equals the distribution of y/l — uK\(0) 
under P . We therefore obtain 



AaK 1 (y) 2 



\u 



u 



E u 



= 4a(l-n)if!(0) 2 



<E U 



(14) 



To estimate the second factor, we desintegrate with respect to the times 
(U\, U2) at which the path of the Brownian bridge hits the position — y for 
the first resp. the last time (on the event that it hits this position at all). 
Again by the strong Markov property and scaling, the distribution K\{—y) 
under P y [.|{i7i = u\,U% = u 2 }] equals the distribution of y/vq — U\ K\ (0) 
under P°. This leads to an estimate analogous to flH}, and allows to conclude 



e 4aifi(j,) 2 " 


< E° 


e 4a*ri(0) 2 ' 


, E v 


e 4aKi(-y) 2 " 


< E° 


^4a^i(0) 2 " 

















By a result due to Levy (see formula (11) in [U]), K\(0) has under 
Raleigh distribution, i.e. 



(15) 



\K 1 (Q)>£) 



_1«2 

e 2 . 



This means that iff (0) is exponentially distributed, and hence, for suitably 
small 5 > 0, E° \ e SKl( - ^ 



is finite. Now (fTTl) follows from (TTHD and (USD. 



So far we have proved existence of a weak solution to © . Weak unique- 
ness is easier to prove, since uniqueness is a local property. Let H be a 
solution to equation ©, and for all n > 1 let T n denote the stopping time 

T n := inf{r > : L r (H r ) > n}. 

By a Girsanov transformation we can change the measure P into a measure 
P under which, for all n G N, the restriction of the process H to the interval 
[0, n A T n ] is standard Brownian motion reflected above 0. Since P and P 
are mutually absolutely continuous, the law of {H sAnA T n , s > 0} under P is 
uniquely determined, for each n > 1. Uniqueness of the law of H solution 
of (El) then follows, since T„ — > 00 a. s. as n — > 00. 
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4.2 Proof of Proposition [3] 

As a by-product of our proof, we will see that the stopping time S x defined 
in (J?]) has finite expectation. A more direct argument for this would make 
the proof of Proposition [3] even shorter, see the discussion in Subsection 15 .11 
Since we have not been able to prove this directly, we circumvent this by 
reflecting the process H z below the level K, and then let K tend to oo. 
To be specific, for K > 0, let H K be the solution of the SDE 

H* = B S + hf (0) - hf (i^), s > 0, (16) 

where L K denotes the local time of H K and B is standard Brownian motion 
defined on the probability space (O, J 7 , P). In other words, H is Brownian 
motion reflected inside the interval [0, if]. 
Let us first note that if we define 

= inf{s > 0, Lf (0) > x}, (17) 

the next result follows readily from Lemma 2.1 in Delmas [2]: 

Lemma 4 For any K > the processes {L Sx (t), 0<t<K} and {L* K (t), 
< t < K} have the same distribution. 

(The intutive explanation of this lemma is as follows: Consider an arbitrary 
level K > 0. The law of Brownian motion reflected in [0, K] equals the law 
of Brownian motion reflected above 0, from which we chop out the pieces of 
trajectories which exceed the level K.) 
We next define the martingale 



M? 



^ 7 {z«)+Lf(^)} 



dB r . 



The same arguments as those in the proof of Proposition [2] show here also 
that for all s > 0, 

Eexp (m?-~{M k ) 8 J =1. 
Therefore there exists a probability measure ¥ K such that for all s > 0, 



dl 



5 A' 



, - oxp(Mf -\{M K ) 
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From Girsanov's theorem, under P , H is a solution of the reflected SDE 

H« = B s + 9 -s- 1 J\z{H?) + L?(H r K )]dr + (0) - h?(K~), s > 0. 

(18) 

We will require 
Lemma 5 

E K [Sff] < oo. 

Proof: We will prove this by a comparison argument. To this end let, 
under P , H K be the solution of 

R* = B s + 6 -s+ l -L« (0) - l -L?(K-), s > 0, (19) 

with L K denoting the local time of R; in other words, R K is a Brownion 
motion with constant upward drift 6/2, reflected above and below K. 
Obviously, for all n = 1, 2, . . . , a 6 [0, if] and e > 0, 

P(Zf +1 (0) - (0) > = a) > P(Zf +1 (0) - I* (0) > e|fl^ = if) := 5, 

where 5 depends on e and K but not on n and a, and is positive at least for 
sufficiently small e. Choosing this e and 5, we see that the probability that 
L (0) increases in the time interval [re, n + 1] by at least e, bounded from 
below by 5, independent of the past. This implies 

E K [Sjf] < oo, (20) 

where S%f is defined by (fT7|) . there with Lf (0) replaced by Lf (0). From a 
classical comparison theorem for SDEs, see e.g. Theorem 3.7, chapter IX of 
[10] , we conclude that < Rg, s > 0, a.s. This implies that 

If (0) < if (0), s > 0, a.s. 

Consequently, 5^ > S 1 ^ a.s.; hence the assertion follows from (j20|) . ■ 
The next subsection will be devoted to the proof of 

Proposition 4 For any K > 0, £/ie process {L% K (t), t > 0} is under F K a 
solution of equation ([2]), killed at time K . 

From this together with Lemma [H Proposition [3] is immediate. 
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4.3 Proof of Proposition [4] 

In this section, x > and K > are fixed. We work under ¥ K and take 
advantage of some of the techniques from [Tj- 

Tanaka's formula yields for any r > and < t < K the identity 
(H? - t)~ = (-*)- - j* ! {H K< t} dH? + ~Lf (t). 

With r := (which is finite P^-a.s. due to Lemma [5]) this yields 

L* K (t) = 2 l {H K< t} dH* . (21) 
Jo 



Plugging (|18p into (|2ip we arrive at 



^fx(t) = x + 2 [ ^ l {H K< t} dB s 
Jo 

+ / l { Hf<t}{0-^{z(Hf) + Lf(Hf)})ds. (22) 
Jo 

It follows from the occupation times formula (see e.g. |10j VI. 1.7) that 

l { H«<t}{0-^z(H?))ds= f\e-2 7 z(u))L^ K (u)du, (23) 
Jo Jo 

while from a generalization of the same formula (see Exercise 1.15 in Chapter 
VI of [10]) we have 

27 I'* l {H K< t} Lf{H«)ds = 2 1 f Lf(u)dLf(u)du 
Jo Jo Jo 

= 7^ [L^{u)) 2 du. (24) 

Let us now abbreviate 

N t :=2 l {H K< t} dB s , 0<t<K. (25) 
Jo 

In order to check that this is a martingale with the appropriate quadratic 
variation, we define, following [7], for all < t < K and s > 



A(s 



J(s,t) r- 



t) : = / 1 {HK< t }dr, r(r,t) := inf{s : A(s,t) > r}, 
Jo 

[ l{H?<t}dB r , £(r,t) := J(T(r,t),t). 
Jo 
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For fixed t, the process £(.,£) as a Brownian motion, arising through a time 
change from the continuous martingale J(.,t). Write F(.,t) for the filtration 
generated by and E% := F(oo,t) for the u-algebra generated by the 

T{s, t), < s < oo. With these slight modifications of the definitions given 
in [7] p. 273, we can carry over all the steps in the proof of [7], Theorem 1, 
to our situation. We will explain here the main ideas and a few details. 

A crucial observation is that every bounded measurable random vari- 
able F can be represented as an ltd integral 



roc roo 

F = K[F] + / v r d£(r,t) = E[F] + / v A{s , t) l {Ha < t} dB s 
Jo Jo 



(26) 



for some F(., i)-predictable process v such that E f Q °° v 2 dr < oo. 

Let u s := 2-l| 0<s<5 xj. This process is predictable, and Lemma[5]implies 
that E Jq 00 u 2 s ds < oo. Moreover for each t > 0, the process 

u(s,t), s > is t)-predictable, 

since u(s,t) = 2 • l{o< s <A(s K ,t)}' an d A(*S^, t) is a £)-stopping time. 
Consequently (recall (|23|) ) 



/'OO /'OO 

N t = J u s l {H K< t }dB s = J u(s,t)d£(s,t) 
is ^-measurable, as well as 

roo f'OO 

° t = Jo u ^ 1 {^^*> ds = J ^ s ^) 2ds - 
Now writing Nt+h ~ Nt as the integral 



u sl{t<H s <t+h}dB s , 

we see from (J5gJ) that E[(7V" m - iV t )F] = 0, in other words, {JV t , < i < iT} 
is an (£t)-martingale. 

Applying Ito's formula to the process s i-> (f Q s u r l{t<H r <t+h}dB r ) , 
s > 0, we obtain as in [7] that 

/•oo /-Sf 

E[(N t+h -N t ) 2 F] = E[ / ^.^ = 4E[ / l {t<Hs < t+h} ds-F], 

Jo Jo 
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which reveals the quadratic variation of (N t ) as (N) t = 4 J Q x 1{#k< 4 } ds. 
Again by the occupation times formula, this equals 2 J* * Lg K (u) du. Conse- 
quently, there exists a Brownian motion {Wt, t > 0} such that 

N t = 2^ y/L^(u)dW u , 0<t<K. (27) 

The proof of Proposition [3] is now completed by combining ([22]) , ([23]) , ([23]) 
P5]l and ([27]). ■ 

4.4 Completion of the proof of the Theorem 

It follows from the description of the law of {Zf , t > Oj^x) made in Section [2] 
that Z is Markov (as a process indexed by x, with values in the set of 
continuous paths from R + into R + with compact support). The fact that 
{Ls x (t), t > 0}{ a .> } enjoys the same property follows from the fact that the 
process Hf := Hs x + r solves the SDE ([6]) with z(t) = Ls x (t) and a Brownian 
motion B which, from the strong Markov property of Brownian motion, is 
independent of {Ls x (t), t > 0}. 

Hence it suffices to prove that for any < x < x + y, the conditional law 
of Ls x+y (-) given Ls x {-) equals that of Z x+V , given Z x . Conditioned upon 
Ls x (-) = z(-), Ls x+y (-) — Ls x (-) is the collection of local times accumulated 
by the solution of ([6]) up to time S y , i. e. it has the law of the process 
{£§ (*), t > 0}, while conditionally upon Z x = z(-), the law of Z x+y - Z x 
is that of Z y,z , solution of equation ([2]). Thus, the assertion of the Theorem 
follows from Proposition [3] 

5 Concluding remarks 

5.1 A possible shortcut in the proof of Proposition [3] 

As a direct consequence of Proposition [3] and the occupation times formula, 
the stopping time S x defined by ([7]) obeys 

s x i r z^dt, 

Jo 

where Z x ' z is the solution of ([2]). This together with a representation of 
J °° Z x '°dt as the random time at which an Ornstein-Uhlenbeck process first 
hits (see [4j) proves 
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Lemma 6 For any x > 0, the stopping time S x defined in (|7|) has finite 
expectation. 

If we could prove Lemma[6]directly from the SDE ([6]), then we could simplify 
our proof of Proposition [3J avoiding the reflection below the arbitrary level 
K. Here is an attempt of a direct intuitive explanation why Lemma [6] holds. 
While climbing up, the Brownian motion with positive drift 9/2 accumulates 
local time at various levels. Sooner or later, it accumulates so much local 
time around some level in M + that the process H governed by © starts to go 
down. It then continues to accumulate local time at various levels, and goes 
back to zero. After reflection at zero, the next excursions will have already 
a stronger drift downwards that awaits H. Remarkably, the recurrence of 
H to the state holds independently of the relative constellations of the 
positive parameters 6 and 7. 

5.2 A second Ray-Knight theorem for Brownian motion with 
a local time drift 

The equation ([6]) is of the form 

H S = B S + h s (0) + I g(H r , L r (H r ))dr, s > 0, (28) 
z Jo 

The proof of Proposition shows that {Ls x (t),t > 0} satisfies the SDE 

Z t = x + [ f(u, Z u )du + 2 f ^fz~ u dW u (29) 
Jo Jo 

with f(t,£) = j g(t,y)dy, provided g is such that (|28|) and ([29]) have unique 
weak solutions which arise via Girsanov transformations from the distribu- 
tions with g = 0, and provided S x = infjs > 0, L s (0) > x} is finite a.s. 
This more general problem will be the object of a forthcoming paper, where 
we will in particular make precise the interaction inside the population, at 
the discrete population level, leading to the continuous limit (|29p . 

Acknowledgement: We thank Jean-Frangois Le Gall for drawing our at- 
tention to the paper [7], Yueyun Hu for helping us streamline the proof of 
Proposition H[ and a referee for a careful reading of a first version that led 
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